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LINEAR  AND  NONLINEAR  THEORY  OF  THE  DOPPLER-SHIFTED 
CYCLOTRON  RESONANCE  MASER  BASED  ON 
TE  AND  TM  WAVEGUIDE  MODES 

I.  Introduction 

During  the  past  twenty  years  the  fast -wave  Cyclotron  Resonance  Maser 
(CRM)  has  been  the  subject  of  considerable  interest  as  an  efficient,  high 
power  source  of  millimeter-wave  radiation.  The  versatility  of  the  interaction 
is  evidenced  by  the  variety  of  device  concepts  which  have  been  developed, 
including  gyrotron  oscillators  and  ampifiers,  gyro-klystrons,  and  high- 
harmonic  gyrotrons.  In  these  devices  the  radiation  frequency  is  close  to  the 
cyclotron  frequency  or  a  harmonic  and  the  beam  energy  is  in  the  weakly 
relativistic  regime.  The  application  of  the  CRM  interaction  to  higher  energy 
beams  and  to  the  generation  of  higher  frequency  radiation  has  led  to  interest 
in  configurations  based  on  the  doppler-shifted  CRM  Interaction.  Potential 
advantages  of  such  configurations  include  a  reduced  magnetic  field 
requirement,  larger  transverse  dimensions  of  the  circuit,  reduced  electron 
beam  pitch  angle,  high  efficiency  of  order  25S,  and  increased  output  power 
from  operation  at  high  voltage. 

The  theory  of  the  CRM  is  quite  well  developed  for  the  weakly  relativistic 
gyrotron  regime  but  the  relativistic  doppler-shifted  regime  has  received  less 
attention.  The  theoretical  approaches  used  most  often  are  the  kinetic  theory 
based  on  the  linearized  Vlasov  equation  and  the  single  particle  theory.  The 
former  leads  to  dispersion  aquations  for  the  Interaction  and  linear  growth 
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rates  or  thresnold  currents.  The  latter  leads  to  both  linear  and  nonlinear 
results  for  the  interaction  parameters.  The  advantage  of  the  kinetic  theory 
is  that  it  can  treat  more  general  beam  distributions  and  usually  involves 
fewer  a  priori  assumptions  than  the  single  particle  method.  In  addition  to 
providing  the  saturated  efficiency,  the  single  particle  theory  can  be  reduced 
to  a  compact  set  of  equations  useful  for  scaling  studies. 

This  paper  presents  a  comprehensive  theory  of  the  Cyclotron  Resonance 
Maser  (CRM)  interaction  in  a  waveguide.  The  circular  waveguide  geometry  was 
chosen  because  it  corresponds  to  the  experimental  geometry  most  often  used  at 
present.  The  kinetic  theory  is  used  to  derive  the  dispersion  relationships 
for  both  TE  and  TM  modes.  The  TE  mode  case  has  been  investigated  by  several 
authors,  but  there  has  been  comparatively  little  work  on  the  TM  mode  case. 
However,  the  TM  mode  interaction  competes  effectively  with  the  TE  mode 
interaction  at  relativistic  electron  energies.  The  conditions  for  maximum 
temporal  and  spatial  growth  rates  are  shown.  The  TM  mode  growth  rates  are 
found  to  vanish  when  the  RF  wave  group  velocity  equals  the  beam  axial  velocity 
("grazing  incidence"). 

The  single  particle  theory  is  used  to  derive  a  compact  set  of  self- 
consistent  nonlinear  equations  for  the  TE  and  TM  mode  interactions.  These 
equations  are  particularly  appropriate  for  the  Cyclotron  Auto-Resonance  Maser 
(CARM)  regime  but  applicability  extends  to  other  regimes  as  well.  The 
conditions  for  optimum  efficiency  are  Investigated  for  oscillator  and 
amplifier  configurations  at  the  fundamental  and  low  order  harmonic 
interactions.  In  the  case  of  a  beam  with  delta  function  distributions  in 
position  and  momentum  the  single  particle  results  in  the  small  signal  limit 
are  shewn  to  be  equivalent  to  the  kinetic  theory  results. 


There  is  an  extensive  literature  on  the  theory  of  Cyclotron  Resonance 
Masers  and  an  exhaustive  review  will  not  be  attempted.  The  initial 
theoretical  work  on  the  radiation  amplification  mechanism  was  presented 
independently  by  Twiss  (1958),  Schneider  (1959)  and  Gaponov  (1959).  A  single 
particle  theory  applicable  to  the  gyrotron  configuration  was  given  in  an  early 
review  by  Gaponov  et  al  (1967).  An  early  application  of  the  kinetic  theory 
was  carried  out  by  Ott  and  Manheimer  (1975)  for  a  planar  configuration.  A 
nonlinear  single  particle  theory  for  the  planar  configuration  was  developed  by 
Sprangle  and  Manheimer  (1975)  and  further  investigated  by  Sprangle  and  Drobot 
(1977).  Application  of  the  kinetic  theory  to  TE  modes  in  cylindrical 
waveguides  was  given  by  Chu  et  al  (1980).  The  TE  and  TM  mode  interactions 
have  been  considered  for  a  beam  with  axicentered  orbits  by  Lau  (1982).  A 
theory  of  ultrarelatlvistic  cyclotron  self-resonance  masers  was  given  by 
Petelin  (1973).  The  single  particle  theory  of  CARM  amplifiers  and  oscillators 
has  been  reviewed  by  Bratman  et  al  (1981,1983).  A  theory  of  doppler  shifted 
relativistic  CRM  amplifiers  has  been  given  by  Ginzburg  et  al  (1981).  Doppler 
shifted  CRM  oscillator  configurations  have  been  investigated  by  FCanavets  and 
Klimov  (1975)  and  by  Vomvoridis  (1982).  The  excitation  of  TM  modes  in  CRM 
oscillators  has  been  investigated  by  Abubakirov  (1983).  An  induced  resonance 
electron  cyclotron  (IREC)  quasi-optical  maser  has  been  investigated  by 
Sprangle  et  al  (1985).  Simulation  studies  of  the  CARM  amplifier  have  been 
carried  out  by  Lin  (1984)  and  by  Lin  and  Lin  (1985). 

The  equations  and  results  presented  in  this  paper  should  provide  a 
useful  starting  point  for  developing  device  concepts  based  on  the  doppler- 
shifted  cyclotron  maser  interaction  and  for  interpreting  on  going 
experiments.  Although  aspects  of  the  present  work  have  been  treated 


previously,  this  work  extends  previous  analyses  by  treating  TE  and  TM  modes  on 
an  equal  footing,  considering  arbitrary  mode  indicies  and  beam  guiding  center 
radius,  and  harmonic  number.  In  addition,  derivations  of  these  theoretical 
results  are  not  generally  available.  Finally,  this  work  provides  a  starting 
point  for  cosidering  more  complicated  systems  such  as  plasma  filled 
waveguides. 

II.  Dispersion  Relationships  from  Kinetic  Theory 

The  present  analysis  considers  an  electron  beam  drifting  in  a  vacuum 
waveguide  with  an  applied  axial  magnetic  field.  Dispersion  relationships  for 
TE  and  TM  modes  in  a  circular  waveguide  are  derived  using  linearized  Vlasov 
theory.  The  derivation  follows  the  approach  of  Chu  and  Dialetis  (1985)  in 
analyzing  the  gyromagnetron.  The  calculation  assumes  the  electron  beam 
current  density  is  sufficiently  low  to  neglect  self-fields.  The  calculation 
is  based  on  the  linearized  Vlasov  equation 

df  (r,p,t)  -  e  (E  +  v  x  B)  .  $♦  f  (r,p,t)  (II. 1) 

P  o 
dt 

where  fQ  is  the  unperturbed  distribution  function  to  be  constructed  from 
constants  of  the  motion  in  an  applied  axial  magnetic  field,  and  E  and  B  are 
the  vacuum  waveguide  rf  fields.  The  transverse  components  are  given  by 


(II. 2b) 


-  Re{T(2)b(?)e~i<ilt} 

,  +  + 

where  e  and  b  are  transverse  vector  functions  given  by 

e  -  2  x  (II. 3a) 


b  ■  - 


(II. 3b) 


for  a  TE  mode,  and  by 


e  «  -  $  i|i 


(II. Aa) 


b  -  -  2  x  ^41 


(11.4b) 


for  a  TM  mode.  The  scalar  function  <|>  satisfies  the  Helmholtz  equation 


♦  kj)*  -  o 

and  the  boundary  condition 


(II. 5) 


I*  -  0 

3n 


(II. 6) 


on  the  waveguide  boundary  for  a  TE  mode,  or  the  boundary  condition 


ijt  •  0 


(II. 7) 


at  the  waveguide  boundary  for  a  TM  mode.  In  Eq.  (II. 5)  denotes  the 
transverse  wavenumber  and  in  Eq. ( II . 6 )  3/3n  denotes  the  normal  derivative. 

5 


The  vector  functions  satisfy  the  following  orthonormality  condition  when 
integrated  over  the  waveguide  cross  section: 


f  ds  h. . 


(II. 8) 


where  h  denotes  either  e  or  b  and  5^  -  1  if  i-j,  6^-  0  if  i  -  j.  The  axial 
field  components  are  given  by 


-  Re{  lk^n(z)  »|>e~iu)t} 

Z  p 

Hi 


(II. 9) 


for  a  TE  mode,  and  by 


2,  2 

c  k 

E  -  Re{  i - -  T(z)1te~lu>t} 

Z  ui 


(11.10) 


for  a  TM  mode.  The  axial  profile  functions  n ( z )  and  T(z)  are  related 
according  to 


T(z) 


i  dll(z) 
oi  dz 


(II.  11  ) 
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-j  hiAcli  i 


for  a  TE  mode  and 


n(z)  -  -  g(,) 

at  dz 


(11.12) 


for  a  TM  mode.  In  the  case  of  a  circular  waveguide,  the  scalar  function  <|>  is 
given  by 


*_<r)  ”  Cmr,Jmn(kmnr)e 
mn  mn  inn  mn 


(11.13) 


where  Jm  is  a  Bessel  function  of  the  first  kind  and  -  xmn/rv’  xmn  is  the 

nth  zero  of  J'  (prime  denotes  differentiation)  for  a  TE  mode  or  the  nth  zero 
m 

of  Jm  for  a  TM  mode,  and  rw  is  the  waveguide  wall  radius.  The  normalization 
constant  is  given  by 


Cmn  ’  {[ir(xL  "  a,2)]1/2jm(xmn)}”1 
am  mn  m  mn 


(11.14) 


for  a  TE  mode  and  by 


C  -  {/if  x  J'(x  ) } 
mn  mn  m  mn 


(11.15) 


for  a  TM  mode. 
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To  evaluate  the  RHS  of  Eq.(II.l),  the  unit  vectors  e  -  p  / [ p  | 

♦  ♦  ♦ 

and  e  »  e  x  e„  are  introduced.  In  terms  of  the  notation  of  Fig.  1 ,  these 
V  U  z 

unit  vectors  satisfy  the  relationships 


cos($-9)e„  *  sin($-9)efl 
r  y 


(II .  1 6a) 


e.  -  -sin($-9)e  +  cosU-9)e. 

<P  r  o 


(II. 16b) 


Using  these  unit  vectors  in  Eq.(II.I)  leads  to 


dt 


P,  3f_ 

e{E  cos(#-9)  ♦  E.sin((*-9)  *  — [B  s in C 4>-9 )  -  B.cos( $-9)  ] }  ^ — 
r  9  •  Yin  r  9  3p. 

o  t 


ep  3f 

[B  oos(H)  -  B  sin(4>-9} 

Ymo  9  r  3pz 


*e{-Esin($-9)  *  EQco3(^-9)  ♦  —  [3Qsin('$-9)  +  B  cos($-9)] 
r  o  Tin  a  r 

o 


-B  f£  ,  1.  !£° 

z  Ymo  p  3<& 


(11.17) 


for  a  TE  mode  and  a  similar  expression  for  a  TM  mode.  Using  Graf's  addition 
theorem  (Abramowitz  and  Stegun  196*0,  and  considering  the  interaction  with  a 
forward  propagating  wave  (H(z)  -  noexp( ikzz) ) ,  the  expressions  in  brackets  in 
Eq. (11.17)  can  be  rewritten  as: 


3 


E  cos($-9)  *  S,sin($-9)  -  He{n  C  k  e  "z 
r  9  o  am  ton 


-Kut-k  z-B0  ♦  »£>  T  j  (k  r  ) 

**  n  mn  r\ 


q  am  o 


J  *  (k  r  )®iq*} 
o+q  mn  L 


(II.l8a) 


-E%in(s-9)  *  Etcos ( $—9 )  -  Refill  C  k  e~iUt“kzz“in0+m2) 
r  9  o  mn  mn 


I  J„(km„rft>Jm+Q(kmnrL)elq0}  <II-l8b) 
k  r.  q  mn  o  m+q  on  L> 

q  on  L  w 


B>3in(<|»-9)  -  B*C03U-9)  -  Re{ - -  noCmnkmne~i  ((*,t"kZZ“me+°2) 

r  8  ca  o  am  mn 


i  <k„o  >eiq*}  (11.13c) 

q  on  o  m+q  mn  u 


3*Sin($-9)  ♦  3*COS(#-9)  -  Se{ - -  ^nCmnkmne  ^  ^ 

9  r  u  o  mn  on 


I  f^r  Jrt(kmnPa)Jm+^kmnrL)  (II*l8d) 

k  r.  q  mn  o  m+q  on  L 

q  mn  L 


k  it 

_  _  mn  _  _  ,  -i(wt-k  z-m9  *  m?) 

3  -  Red -  II  C,  k  e  z  2 

z  to  o  mn  mn 


l  J(i(kmnrn)  JB+fl(ktnnr’L)  e  ' 

^  q  mn  o  m+q  mn  L 


(II.l8e) 


The  unperturbed  distribution  function  fQ  is  now  treated  as  a  function  of  rQ, 


and  ?>  and  the  partial  derivatives  are  written  as 


(■>11  WTym’WaV'V 


-  .  v%  1%  .  - 


£  .  £  .  J_<.M  .  .-15,  £ 

3pt  3pt  Ve  3ro 


(II. 19a) 


1  3fo 


pt30  2moQe 


-  .'li)  £ 


(11.19b) 


where  In  the  above  equations  tQ  -  fQ(pz,pt,*)  on  the  LHS  and  fQ  »  fQ(pz,pt ,rQ) 


on  the  RHS.  Substituting  Eqs.(II.l8)  and  (II. 19)  into  Eq. (11.17)  yields 


Z  “  Re{f%  C  k  l  e 

~r~  -  mn  mn 

at  ui  s 


-  i  ( ait-kzz-m$-  ( s-tn )  i+o|) 


kzpz,  sfo  .  Vt  5fo. 


“.-.'Vo1  C<“  •  ’fif1  157  *  TT  W) 

O  u 


*o  5pz' 


-  -4-  (j*  (k  r  )J  (k  r,  )  (•> 
in  8  3-m  on  o  s  on  L 


~  -kJ  <kmOJ  (kfn„r’T)  '  J  (kmnrr>)J  (kmnrr)3 

2  _  , on  o  .on  L  _  _ . , on  o  ..on  L 


s-o- 1 


3-m+ T 


kranpt  i^o-i  , 
Yn0  f  3ro 


(11.20) 


In  what  follows,  the  the  last  terra  on  the  RHS  of  Eq. (11.20)  involving 


3fQ/3ro  will  be  neglected.  It  can  be  shown  that  for  configurations  which 


optimize  the  Cyclotron  Maser  interaction,  contributions  arising  from  this  term 


are  negligible.  Instead  this  terra  leads  to  the  Peniotron  interaction  (Dohler 


et  al  (1978),  Vitello  (1984),  Zhang  (1985)).  Eq. (11.20)  is  easily  solved 


using  the  method  of  characteristics  (Krail  and  Trivelpiece  1973)  in  which  the 


RHS  is  integrated  with  respect  to  time  along  the  unperturbed  orbits, i.e., 


z(t')  -  z  ♦  v(t'-t),  $(t*)  -  *  *  B  (t'-t)/Y,  H(t')  -  ♦  ♦  8  (t '-t)/Y. 
zo  e  e 


For  a  TE  mode  this  gives 
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Re{e  ~  C  k  I  fIzCp  ,pt ,p  ) 
oj  mn  mn  am  o  z  z 
a 


.  -i(wt-  k  z-m$-(s-m)$+mr) 
ie  z  2 


«-k  v  -  sfl  /Y 
z  z  e 


(11.21 ) 


where 


J  )C(“ 

mn  o  s  mn  L 

3-0 


Vz,  %to  .  w  3f< 


2  2'  _ O  +  2  t  _ 0-j 

Ym0  3pt  Ym0  3p2J 


(11.22) 


and  for  a  TM  mode 


"•‘-"VwAm  ?  WW 


-i(<ut-k  z-m$-(3-m)$+tfc) 

e  z  2 


z  z  e 


(11.23) 


where 


FOT<rO,Pt’pZ) 


ck  3f 
z .  o 


k  r.  Js-m<kinnro,'Vkiiin,"L)  i  c  "  <a  1  3p„ 


f-  J _  .  "anpt  ,  ffo. 

^  Ym  c  afl  me'  3p, 
o  e  o  z 


(11.24) 


In  Eq. (11.24)  -  k  c  and  Q  is  the  non-relativistic  cyclotron  frequency, 

mn  mn  e 


The  dispersion  relationships  are  obtained  from  Maxwell's  Equations.  The 


wave  equation  for  a  TE  mode  in  cylindrical  symmetry  is 


,  1  3  -2'_  r1  3  .  .Oh 

‘7  “2  '  7  ,Bz  •  “<CF  3r  rJ9  > 

C  at 


a,<'> 
lit  j 

r  38  J 


(11.25) 


where  the  components  of  the  perturbed  current  density  are  given  by: 
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Similarly,  the  wave  equation  for  a  TM  mode  is 
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where  the  perturbed  axial  current  and  charge  densities  are  given  by: 
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Substituting  Eq.(II.9)  into  Eq.(II.25),  multiplying  by  rJm(kmnr)e  1  and 
integrating  with  respect  to  r  and  9  leads  to: 
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Integrating  the  RHS  of  £q. (11.29)  by  parts  and  applying  Graf's  addition 
theorem  yields 
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Finally,  using  the  relationship 
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the  dispersion  relationship  for  the  TE  mode  Interaction  with  a  given  harmonic 


s  is  obtained: 
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By  a  similar  analysis  the  dispersion  relationship  for  a  TM  mode  is  obtained: 
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Writing  the  unperturbed  distribution  function  in  the  form 
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(11.34) 
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where  a  is  the  number  of  beam  electrons  per  unit  volume,  and  integrating  the 
o 

RHS  of  Eqs. (11.32)  and  (11.33)  by  parts  with  respect  to  p  and  p  ,  the  TE 

t  Z 

mode  equation  becomes 
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For  a  TE  mode  the  coefficient  is  defined  as 
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The  TM  equation  becomes 
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For  a  TM  mode  the  coefficient  i3  defined  as 
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Sqa. (11.35)  and  (11.38)  can  be  simplified  in  the  case  of  a  cold,  thin 
annular  beam.  This  corresponds  to  the  delta  function  distribution  function 
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where  is  the  area  of  the  beam  cross  section  and  r&  is  the  average  beam 
radius.  It  is  also  convenient  to  introduce  the  plasma  frequency 
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and  the  beam-waveguide  fill  factor 


irr 

w 


which  is  usually  much  less  than  unity.  Substituting  Eqs. (II. 41  Mil. 43)  into 
Eqs.(II.35)  and  (II.38)  leads  to  the  cold  beam  dispersion  equations  for  the  TE 
mode: 
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and  for  the  TM  mode: 
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Eq. (11.44)  agrees  with  the  result  derived  previously  by  Chu  et  al  (i960) 
except  that  their  result  includes  additional  terms  in  the  expression  for 


Qsn.  These  terms  are  not  included  in  the  present  analysis  because  the  term 


involving  3fQ/3ro  in  Eq. (11.20)  has  been  neglected.  In  conf igurations 


optimized  for  the  cylcotron  maser  interaction  these  terms  are  small.  However 
in  a  related  device  called  the  peniotron  which  involves  large  electron  guiding 


center  drifts  these  terms  can  be  important  (Zhang  1985).  Thus  the  present 
analysis  does  not  apply  to  this  device.  As  discussed  by  Chu  et  al  (1980),  the 
first  term  on  the  RHS  of  Eq. (11.44)  is  the  source  of  the  instability,  while 
the  second  is  a  stabilizing  term.  For  relativistic  beams  with  signifigant 
transverse  velocity,  the  second  term  is  much  smaller  than  the  first  and  can  be 
neglected. 

Comparing  Eqs. (11.44)  and  (11.45)  shows  that  the  TE  and  TM  mode 

dispersion  relations  are  similar  in  form  but  there  are  significant 

differences.  The  TM  mode  instability  tends  to  be  reduced  compared  to  the  TE 

mode  case  by  two  effects:  First,  the  RHS  of  Eq. (11.45)  involves  an  additional 

2 

stabilizing  term  proportional  to  (w  -  k  v  -  sQ  /Y)  .  For  frequencies  near 

Z  ZO  6 

the  doppler-shifted  cyclotron  frequency  this  term  will  be  usuually  be  small 

compared  to  the  others  and  can  be  neglected.  Secondly,  the  other  terms  on  the 

-1  2 

RHS  of  Eq. (11.45)  include  the  factor  (8„  -  8  .  )  which  vanishes  when  the  beam 

z  ph 

2 

velocity  equals  the  wave  group  velocity  since  v  •  c  /vph  .  Thus,  the  TM  mode 
instability  vanishes  when  the  vacuum  waveguide  mode  and  beam  mode  are  at 
grazing  incidence.  This  effect  has  been  pointed  out  by  Abubakirov  (1983). 


The  temporal  growth  rates  for  the  instability  at  resonance  are  obtained 
by  substituting  m  -  u  *  Aw,  k  -  k  into  the  dispersion  equations,  where 

O  Z  ZO 

(«  ,  k  )  is  the  point  of  intersection  between  the  uncoupled  waveguide  mode 
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and  the  beam  mode 
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Neglecting  small  terms  as  discussed  above  leads  to  the  cubic  equations 
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for  the  the  TE  and  TM  modes,  respectively.  Setting  Au>  -  Aid  «•  iAd). , 

r  i 

Eqs. (11.48)  and  (11.49)  have  the  solutions 
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respectively,  where  u)  is  the  width  of  the  beam-wave  interaction.  The 

r 

temporal  growth  rate  is  given  by 


Ad),  *  /3  Ad) 


(11.52) 


Inspection  of  Eqs. (11. 50) -{11. 52)  shows  that  for  a  given  electron  beam  and 
waveguide,  the  growth  rate  is  highest  for  a  wave  near  cutoff  and  approaches 


zero  as  the  phase  velocity  approaches  the  free  space  value  for  both  TE  and  TM 
modes.  In  addition,  the  growth  rate  vanishes  at  grazing  incidence  for  TM 
modes  (vz  -  vg). 

As  pointed  out  by  Bratman  et  al  (1983),  this  effect  can  be  explained  by 
considering  a  rest  frame  moving  with  the  group  velocity  of  the  wave.  In  this 
frame  the  mode  is  at  cut-off  (k2  -  o)  and  in  the  case  of  a  TM  mode  has  only  a 
longitudinal  electric  field  component.  If  the  electron  beam  drift  velocity 
vz0-vgr,  then  in  this  frame  the  electron  motion  is  purely  transverse  and  the 
rf  field  does  no  work  on  the  electrons.  This  effect  ha3  also  been  pointed 
out  by  Lau  (1982)  for  a  beam  with  axi-centered  electron  orbits. 


The  spatial  growth  rate  can  be  calculated  by  substituting 

u»  -  u»  and  k  -  k  *  Ak  into  Eqs. (11.44)  and  (11.45).  The  result  can  be 
o  o 

expressed  as 
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where  3  »  8  .  .  From  Eqs. ( II . 50 )— (11.53)  the  small  signal  spatial  growth 
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rates  expressed  in  units  of  d3/m  are 
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for  a  TE  mode,  and 
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for  a  TM  mode.  Eqs. (11.5*0  and  (11.55)  do  not  apply  when  the  instability 
becomes  absolute  such  as  occurs  at  the  waveguide  cutoff  frequency.  The  onset 
of  absolute  instability  for  a  TE  mode  at  the  grazing  incidence  condition  has 
been  analyzed  by  Lau  et  al  (1981). 


III.  Single-Particle  Theory 


In  this  section  the  Cyclotron  Maser  interaction  is  treated  in  the  single 
particle  approximation.  A  thin  annular  beam  propagating  in  a  circular 
waveguide  in  the  presence  of  an  applied  axial  magnetic  field  is  considered. 

The  electrons  follow  helical  trajectories  in  the  applied  magnetic  field  about 
guiding  centers  located  at  a  radius  rQ  from  the  symmetry  axis.  Slow-time- 
scale  nonlinear  equations  are  obtained  for  the  interaction  with  a  given 
harmonic. 


The  starting  point  is  the  Lorentz  force  equation  for  an  electron 
drifting  in  an  applied  predominantly  axial  magnetic  field  and  experiencing 
perturbing  rf  fields: 
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where  Bq  is  the  applied  magnetic  field 


( III.  1 ) 
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f  and  B  are  the  rf  fields,  and  f  is  the  relativistic  factor 
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(III. 3) 


To  obtain  slow-time-scale  variables,  the  transverse  momentum  is  expressed  in 
the  form 
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(III. 4) 


where  pt  and  $  are  the  slowly  varying  magnitude  and  phase,  Q  is  the  reference 
cyclotron  frequency 
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(III. 5) 


where  zQ  -  z(tQ ),  and  t  *  t  -  tQ  where  t0  is  the  time  the  electron  enters  the 
interaction  region.  Substituting  Eq.(III.4)  into  Eq.(III.l)  and  neglecting 
drifts  of  the  electron  guiding  center  due  to  the  rf  fields,  the  equations  of 
motion  can  be  expressed  in  the  form 
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In  the  equations  which  follow  the  magnetic  field  gradient  terms  will  be 
dropped  corresponding  to  the  case  of  a  constant  applied  axial  magnetic 
field.  Substituting  Eqs.  (II.2)-(II. 4) ,  and  (II.9)-(II.10)  into  Eqs.(III.6) 
leads  the  following  equations  of  motion  for  a  TE  mode: 
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and  for  a  TM  mode: 
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where 
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3y  using  the  expression  for  the  scalar  function  t| t  given  in  Eq.(II.13) 
in  Sqs.(III.7)  and  (III. 8),  applying  Graf's  Addition  Theorem  for  Bessel 
functions,  and  considering  the  interaction  with  the  sth  harmonic  of  the 
cyclotron  frequency,  the  cyclotron  maser  equations  of  motion  for  a 
circular  waveguide  can  be  written  in  the  form: 
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for  a  TS  mode,  and 
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for  a  TM  mode;  where  the  phase  variable  has  been  replaced  by 
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for  a  TM  mode.  For  a  circular  waveguide  and  a  single  harmonic,  these 
equations  reduce  to 
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respectively. 


As  pointed  out  by  Bratman  et  al  (1981),  the  relationship 
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holds  in  the  case  of  the  Cyclotron  Maser  interaction  for  a  constant 
amplitude  free  space  wave  of  the  form 
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where  1I0  is  a  constant.  By  substituting  Eq. (11.11)  into  (III. 10c), 
using  Eq.(III.l9),  and  comparing  with  Eq. ( III .16),  Eq.(III.l8)  is 
readily  verified  for  TE  modes  in  a  waveguide.  Eq . C III .18)  is  not  in 
general  true  for  TM  waveguide  modes  as  can  be  seen  by  subtracting  vph 
times  Eq. (III. 11c)  from  Eq.(III.l7)  to  obtain: 


[“  -  ~  -  6  ]Re{iT  e"i(A"kz2)}  (III. 20) 
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where  is  the  kinematic  phase  shift  parameter 
.  .  fzo  sfl 

o  '  3  .  “  oi  (III. 21 ) 

ph 

is  the  constant  wave  amplitude,  and  8  -  v  . /c,  B  •  v  /c. 

°  ph  ph  z  z 

Sq . ( III .19)  shows  that  Eq.(III.l8)  holds  for  the  TM  mode  interaction 
for  a  constant  amplitude  wave  only  when  vph  -  c,  the  free  space  limit. 


It  is  convenient  to  introduce  the  following  normalized  momenta  and  energy: 


the  normalized  axial  coordinate 


Z  ■  (tfZ/  c 


(III. 23) 


and  to  introduce  the  slowly  varying  momentum  phase 
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(III. 24) 


Transforming  to  spatial  derivatives,  the  TE  mode  energy  and  phase 
equations  for  a  constant  amplitude  wave  can  then  be  written  as: 
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Using  conservation  of  energy 
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and  Eq.(III.l8),  the  momentum  amplitudes  can  be  expressed  in  terms  of 
the  normalized  energy  as  follows: 
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For  Interactions  with  non-constant  amplitude  waves,  such  as  occurs  in  a 
traveling  wave  amplifier,  Eq.(III.l8)  does  not  hold  for  either  TE  or  TM 
modes.  In  the  case  of  a  TE  mode,  a  slow  axial  variation  of  the  wave 
amplitude  affects  the  axial  momentum  according  to: 


3J(k  r,  ) 
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where  q'  is  defined  by 


q'  •  pj  *  a/V 


(III .29) 


For  a  slow  variation  in  wave  amplitude,  it  is  sufficient  to  include 
the  axial  momentum  correction  q*  in  the  phase  equation  according  to: 
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A  similar  analysis  leads  to  somewhat  more  complicated  equations  for  the 
TM  mode  interaction: 
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In  the  CARM  regime  (Bratman  et  al  1981)  the  phase  velocity  is  close 
to  the  speed  of  light  and  the  following  conditions  apply: 


1  -  i'J.  <<  1 
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Under  these  conditions  Eq.(III.l8)  holds  to  a  good  approximation  for  TM 

modes  so  that  that  momentum  amplitudes  can  be  calculated  from 

Eqs.dII.27) .  Using  this  fact  and  neglecting  terms  proportional  to 
-2 

(1  -  ) / C 1  *  Sj/Spft)*  the  equations  of  motion  for  the  CARM  regime 

become 
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for  a  TE  mode  and 
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for  a  TM  mode,  where  the  effect  of  non-constant  wave  amplitude  is 
neglected  except  in  the  phase  equation.  In  Eq. (III. 33b) ,  q'  is  given 
by  Eq.(III.28).  To  further  simplify  the  equations,  the  following 
variables  defined  by  Bratman  et  al  (1981)  are  introduced: 
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In  this  notation  the  momentum  amplitudes  in  the  CARM  regime  are  given 


p*  -  3t0/l  -  u 


(III. 36a) 
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The  parameter  b  characterizes  how  strongly  the  axial  momentum  and  velocity 
change  with  change  in  electron  energy. 

In  Cyclotron  Maser  interactions  with  the  fundamental  or  low  order 
harmonics  the  electron  Larmor  radius  is  usually  small  compared  to  the 
transverse  dimensions  of  the  waveguide  mode.  In  this  case  the  Bessel 
function  of  order  s  -  1  occurring  in  the  equations  of  motion  may  be 
replaced  by  the  leading  term  of  the  small  argument  expansion  with 
little  loss  of  accuracy.  In  terms  of  the  present  notation 
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for  a  TM  mode.  Using  Eq3.(III.35)  -  (III. 39),  the  CARM  energy  and 
phase  equations  of  motion  may  be  reduced  to  the  following  form: 
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for  TE  modes,  and  to 
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for  TM  modes.  For  non-constant  wave  amplitudes,  the  axial  momentum 
correction  is  given  by 
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for  TE  modes  and  by 
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for  TM  modes. 


The  above  equations  describe  the  motion  of  an  electron  subject  to 


prescribed  rf  fields.  In  a  self-consistent  formulation  the  wave 


amplitudes  should  be  calculated  from  the  induced  AC  current  and  charge 


densities.  Assuming  a  time  dependence  of  the  form  exp(-iwt),  the 


following  wave  equations  are  obtained  from  Maxwell’s  Equations 


appropriate  for  TE  and  TM  modes,  respectively: 
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where  in  Sqs.(III.44)  E  and  E  are  complex  phasors  and  J.  ,  J  ,  and  p  are 

Z  uU)  ZU)  (a) 


the  fundamental  Fourier  components  of  the  AC  current  and  charge  densities. 


Using  the  prescription  outlined  by  Flyagin  et  al  (1977)  and  by  Fliflet  et  al 


(1982),  the  following  wave  equations  for  the  z-dependent  mode  amplitudes  H 


and  T  can  be  derived: 
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where  the  Eq.(45a)  applies  to  TE  modes  and  Eq. (III. 45b)  applies  to  TM 

modes.  Considering  the  interaction  with  a  co-propagating  wave,  i.e., 

II  -  H  exp(ik  z),  T  -  T  exp(ik  z),  neglecting  second  derivatives  of 
O  '  z  o  z 

HQ  or  Tq,  introducing  normalized  variables,  and  approximating  the  Bessel 
functions  of  order  s,  leads  to  the  following  equations  for  the 
normalized  TE  wave  amplitude: 
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and  for  the  TM  wave  amplitude: 
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where  the  normalized  current  parameter  is  given  by: 
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for  a  TE  mode  and  by: 


for  a  TM  mode.  Eqs. (III. 40) ,  (III. 42)  and  (III. 46)  form  a  self- 

consistent  set  for  the  TE  mode  traveling  wave  interaction.  Similarly, 

Eqs. (III. 41 ) ,  (III. 43)  and  (III. 47)  form  a  self-consistent  set  for  the 

TM  mode  traveling  wave  interaction  in  the  CARM  regime.  Although  the 

derivation  has  been  carried  out  for  the  CARM  regime,  it  is  readily  shown  that 

the  above  equations  apply  in  the  gyrotron  limit  as  well  (b  •+  o,  8  .-*■  «). 

ph 

It  is  of  interest  to  compare  the  results  of  single  particle  theory 
with  the  results  of  the  kinetic  theory  derived  above  for  a  tenuous  beam 
with  5  -  function  distributions.  Expressing  the  wave  amplitude  in  the 
form 
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(III. 50) 


solving  the  single  particle  equations  in  the  small  signal  limit  and 
averaging  with  respect  to  the  initial  electron  phase  leads  to: 
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When  the  beam  current  parameter  is  small,  Eqs.(III.5l)  and  (III. 52)  reduce  to 


(T  -  A)  r2  -  -  |  (III. 5*) 

which  -  as  pointed  out  by  Bratman  et  al  (1981)  -  has  the  same  form  as 
the  well-known  Traveling- Wave  Tube  amplifier  dispersion  equation 
(Gewartowski  and  Watson,  1965).  The  normalized  growth  rate  at 
resonance  (  A  -  0)  is  then  given  by 

/r  i  1/3 

Imr  -  -f-  [-J-]  (III. 55) 

for  either  TE  or  TM  modes. 

By  substituting  unnormalized  parameters  into  Eqs.(III.55)  and 
(II 1.62)  and  using  the  small  argument  approximation  for  the  appropriate 
Bessel  functions  in  Eqs. (11.50)  -  (11.51),  it  is  readily  shown  that  the 
3ingle  particle  TE  and  TM  mode  growth  rates  agree  with  the  kinetic 
theory  results  given  in  Eqs. (11.50)  -  (11.53). 

In  terms  of  the  present  normalized  parameters,  the  electronic 
efficiency  is  given  by: 

8 2  <u(y)> 

n  -  - -  (111.56) 


where  u  is  the  point  where  the  interaction  terminates  and  <  >  denotes 

averaging  with  respect  to  the  initial  momentum  phase  angle.  The  total  time- 
averaged  power  of  a  traveling-wave  integrated  over  the  waveguide  cross  section 
at  any  point  is  given  by: 
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In  a  high  Q  oscillator,  the  RF  field  axial  profile  is  essentially 
determined  by  the  cavity  and  not  by  the  interaction  with  the  beam.  In 
this  case  the  equation  for  the  field  amplitude  can  be  neglected  and  the 
efficiency  can  be  found  by  integrating  the  equations  of  motion  in  a 
prescribed  cavity  RF  field.  The  required  beam  power  is  found  by 
applying  the  power  balance  equation 
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where  VQ  is  the  beam  voltage,  Q  is  the  cavity  Q  factor  neglecting  wall 
losses,  and  W  is  the  cavity  stored  energy  which  is  given  by 
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For  a  cavity  defined  by  a  uniform  waveguide  terminated  by  sections  of 
high  reflectivity,  the  field  amplitude  is  uniform  and  the  stored  energy 
is  given  by: 
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The  snail  signal  efficiency  for  a  uniform  field  amplitude  and  the  scn 
harmonic  is  given  by: 
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and  L  is  the  interaction  length.  The  parameter  corresponds  to  the 
parameter  A  defined  by  Chu  0978).  Considering  the  small  signal 
limit,  substituting  Eq.(III.6l)  into  Eq.(III.58),  and  transforming  to 
normalized  variables  leads  to  the  following  expression  for  the 
threshold  beam  current  for  either  TE  or  TM  modes: 
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The  Q  factor  for  a  cavity  with  unifora  wave  amplitude  is  given  by 
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where  R1  and  R2  are  the  reflection  coefficients  at  the  cavity  input  and 
output.  Substituting  Eq.(III.64)  into  (III. 63)  leads  to 
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IV .  CALCULATIONS 

A.  High  Q  oscillator 

Applications  of  the  present  theory  have  been  considered  by  several 
authors,  therefore  only  selected  calculations  are  presented  here. 

In  the  case  of  a  high-Q  oscillator  with  a  uniform  axial  magnetic 
field,  the  electrons  interact  with  an  approximately  constant  amplitude  wave 

A 

and  the  normalized  efficiency  n  -  <u>  depends  on  only  five  parameters:  the 
wave  amplitude  FQt  the  interaction  length  u,  the  detuning  parameter  A,  the 
harmonic  number  s,  and  the  parameter  b  which  characterizes  how  strongly  the 
electron  longitudinal  velocity  changes  with  a  change  in  electron  energy  during 
the  interaction.  A  theory  of  the  quasi-optical  gyrotron  which  corresponds  to 
b  -  0,  involving  four  parameters  was  obtained  by  Bondeson,  Manheimer  and  Ott 
(1983).  In  the  following  calculations  the  detuning  parameter  is  treated  as  an 

A 

optimization  variable.  Figure  2  shows  the  normalized  efficiency  n  as  a 
function  of  b  for  several  values  of  F0.  The  calculations  are  for  s-1  and  are 
optimized  with  respect  to  y.  The  corresponding  values  of  u  are  given  in 


u  tu, 


Figure  3*  Figure  4  compares  efficiency  of  the  first  four  harmonics  for  FQ  - 
0.2  and  optimized  A  and  u.  The  optimum  interaction  lengths  for  the  higher 
order  harmonics  (not  shown)  are  only  slightly  different  from  the  results  for 
s»i .  The  actual  efficiency  is  obtained  for  these  results  by  using  Eq. 

(III. 56).  Note  that  the  results  given  in  Figures  2  and  3  are  independent  of 
mode,  electron  energy  and  orbit  pitch  angle,  and  the  wave  frequency  and  phase 
velocity. 


The  limit  b  -  0  corresponds  to  the  gyrotron  interaction  where  the 
wave  is  close  to  cutoff  and  the  wave  frequency  is  approximately  equal  to  the 
relativisitc  cyclotron  frequency.  This  limit  gives  the  highest  normalized 

A 

efficiency  (  n  -  425)  and  has  been  extensively  investigated.  As  is  well  known 
the  efficiency  for  this  limit  can  be  enhanced  considerably  by  contouring  the 
axial  profile  of  the  RF  field.  Normalized  efficiencies  of  order  80)1  have  been 
reported  (Gaponov  et  al  1981). 


As  discussed  by  Bratman  et  al  (1981,1983),  the  Cyclotron  Auto- 
Resonance  Maser  (CARM)  regime  corresponds  to  b  -  0.5.  The  single-particle 
efficiency  is  optimized  (n  5  0.5)  by  choosing  8  »  1/Y.  A  large  doppler 

SP  w 

2 

upshift  -  ai  S  Y  Q  where  Q  is  the  relativistic  cyclotron  frequency  -  and  good 

C  0 

A 

bunching  efficiency  (characterized  by  n  )  occur  when  in  addition 
-2  -2 

( 1  —  B ph )  <  Y  •  Figure  2  shows  that  a  maximum  normalized  efficiency  of  365  is 
obtained  at  b  -  0.5  for  FQ  -  0.2,  A  -  0.6,  and  y  *  8.  This  result  was  given 
previously  by  Bratman  et  al  (1981).  One  expects  that  considerably  higher 
efficiencies  should  be  achievable  by  enhancement  techniques  similar  to  those 
used  in  gyrotrons. 
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Another  regime  of  doppler-shifted  operation  has  been  investigated  by 


Vomvoridis  (1982)  for  the  oscillator  and  by  Ginzburg,  Zarnitsyna,  and 

Nusinovich  0  981).  In  this  case  3,  -  (1-Y 1  )/B  .  and 

zo  o  ph 

sf  *23  ?Y  ^  ( 1  -Y  1 )  +  (1-3  2)  (1-Y  2).  An  interesting  feature  of  this  choice 
to  ph  o  o  ph  o 

is  both  and  3  can  vanish  simultaneously,  that  is,  the  total  kinetic  energy 
is  available  to  the  interaction  (Vomvoridis  1982)  and  the  single-particle 
efficiency  approaches  100?.  In  this  regime  b  -  1  and,  as  shown  in  Figure  2, 
the  maximum  normalized  (and  unnormalized)  efficiency  is  22?.  As  in  the  other 
regimes,  higher  efficiencies  should  be  achievable  through  enhancement 
techniques. 

Based  on  these  calculations  design  parameters  are  given  in  Table  I 
for  a  100  GHz  CARM  oscillator  experiment  using  600  kV  beam  and  choosing 


B.  Traveling-Wave  Amplifier 

Calculations  have  been  carried  out  to  investigate  a  CARM  regime 
amplifier.  Figure  5a  shows  the  normalized  efficiency  as  a  function  of 
normalized  current  for  b  -  0.454  and  FQ  -  0.0018.  The  results  have  been 
optimized  with  respect  to  A  and  interaction  length  u.  As  pointed  out  by 
Ginzburg,  Zarnitsyna  and  Nusinovich  (1981)  the  saturated  efficiency  and  wave 
amplitude  are  independent  of  the  input  amplitude  when  the  input  amplitude  is 
sufficiently  small.  The  results  shown  in  Figure  5  represent  this  limit.  The 
values  of  saturation  wave  amplitude  and  interaction  length  corresponding  to 
Figure  5a  are  shown  in  Figures  5b  and  5c,  respectively.  Except  at  low  beam 

A 

current  (IQ  <  0.05)  which  re-  ires  smaller  values,  the  optimum  ietuning 
is  A  •  0.4. 
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The  effect  of  increasing  the  input  wave  amplitude  for  a  given  beam 

A 

current  is  shown  in  Figure  6  for  a  normalized  current  of  I  -  0.60  and  b  - 

o 

0.47.  The  results  have  been  optimized  with  respect  to  A  and  u  as  in  Figure 
5.  Figure  6a  gives  the  normalized  efficiency  dependence  on  the  input  wave 
amplitude.  As  noted  above,  the  efficiency  is  essentially  constant  over  a  wide 
range  of  input  amplitude,  then  rises  and  saturates  at  a  power  gain  of  about  10 
dB.  The  values  of  saturation  wave  amplitude  and  interaction  length  are  shown 
in  Figures  6b  and  6c.  The  optimum  detuning  is  A  »  0.4. 


The  design  of  a  high  power  CARM  amplifier  operating  in  the  TE1 1 

circular  waveguide  mode  with  parameters  similar  to  the  NRL  VEBA  free  electron 

laser  experiment  (Gold  et  al  1984}  has  been  considered.  The  design  parameters 

are  given  in  Table  II.  The  beam  energy  is  1  MeV,  the  current  is  500  Amps,  the 

waveguide  radius  is  5.4  mm,  and  the  design  operating  frequency  is  94  GHz. 

Figure  7  compares  the  small  signal  spatial  growth  rates  for  the  TE^  and  TMQ1 

modes  as  a  function  of  wave  phase  velocity  (and  frequency)  in  the  design 

waveguide  (rw  -  0.54).  The  growth  rates  were  calculated  using  the  parameters 

of  Table  II  and  Eqs. (11.54)  and  (11.55)  with  small  argument  expansions  for  the 

Bessel  functions  with  argument  k^r^.  The  guiding  center  radius  was  taken  to 

be  zero  for  the  TE^  mode  and  as  2.1  mm  for  the  TMqi  mode.  Figure  7  shows  as 

discussed  above  that  the  TM01  growth  rate  vanishes  at  the  grazing  incidence 

condition  (3  3  »  0.88).  Except  for  this  effect  for  TM  modes,  the  growth 

ph  zo 

rates  increase  monotonically  as  the  wave  approaches  the  cut-off  condition. 

This  calculation  assumes  the  instability  is  convective,  at  some  point  it 
becomes  absolute  and  the  calculation  is  invalid  (see  Lau  et  al  1981). 
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The  peak  output  power  of  the  design  is  110  MW  at  an  electronic 
efficiency  of  22%  (<u>  -  0.34).  The  normalized  interaction  length  is  u  -  15.4 
for  an  input  power  of  1  kW  (FQ  -  0.0019)  and  the  total  gain  is  50  dB.  The 
small  signal  and  nonlinear  efficiencies  are  plotted  as  a  function  of  the 
detuning  parameter  A  in  Figure  8.  The  small  signal  result  has  been  scaled  by 
a  factor  of  10  for  clarity.  Both  curves  were  calculated  numerically,  the 
small  signal  results  correspond  to  an  input  power  of  1  W  and  the  nonlinear 
results  correspond  to  1  kW.  The  FWHM  linear  and  nonlinear  bandwidths  are  10% 
and  3%.  respectively.  Note  that  the  bandwidth  cannot  be  directly  obtained 
from  Figure  8  because  A  is  not  proportional  to  frequency.  The  bandwidth 
could  be  increased  by  decreasing  the  total  gain.  Figure  8  shows  that  the 
region  of  maximum  nonlinear  gain  is  signifigantly  detuned  from  the  region  of 
maximum  small  signal  gain. 

As  a  final  example  a  2  THz  second  harmonic  (s-2)  CARM  amplifier  i3 
considered.  Operation  at  high  power  at  this  frequency  requires  that  the 
waveguide  transverse  dimension  be  much  larger  than  the  radiation  wavelength 
(A-0.15  mm).  This  can  be  achieved  partly  by  operating  in  a  higher  order  mode 
but  primarily  by  locating  the  beam  line  -  waveguide  mode  intersection  point 
far  above  the  waveguide  cut-off.  This  ca  e  represents  the  "auto-resonance" 
limit,  however,  as  discussed  above  the  rate  of  gain  is  reduced  leading  to  a 
longer  interaction  length  and  increased  sensitivity  to  beam  velocity  spread. 

In  addition,  there  is  an  increased  probability  for  exciting  spurious  modes. 
Nevertheless,  this  design  regime  is  of  interest  due  to  the  lack  of  high  power 
sources  at  THz  frequencies. 


The  waveguide  radius  is  chosen  to  be  0.5  cm  which  allows  an  annular 
electron  beam  of  radius  0.4  cm,  and  the  operating  mode  is  chosen  to  be  the 
TE10  1  whispering  gallery  mode.  This  corresponds  to  a  wave  frequency  to  cut¬ 
off  frequency  ratio  of  w/ui  -  17.8  or  a  wave  phase  velocity  of  1.0016c.  For  a 

CO 

1.2  MeV  beam  interacting  at  the  second  harmonic  and  with  StQ-  the  doppler 

upshift  factor  is  u/sfl  -  10.2.  The  required  magnetic  field  is  115  kG  which  is 
within  the  state-of-the-art  for  superconducting  solenoids.  The  design  was 
optimized  for  a  beam  current  of  1  kA  and  assumes  an  input  power  of  1  kW.  The 
maximum  efficiency  and  power  were  calculated  to  be  19$  and  223  MW  at  an 
interaction  length  of  134  cm.  The  nonlinear  instantaneous  bandwidth  is 
approximately  0.04$  and  the  total  gain  is  53  dB  corresponding  to  an  average 
gain  rate  of  0.4  dB/cm.  The  design  parameters  are  summarized  in  Table  III. 

The  normalized  parameters  for  this  design  are  b-0.518,  Iq-0.0044,  Fq-0.0001 , 
Fgaj. -0.050,  u-92,  A-0.15,  and  <u>«0.28. 

V.  Conclusions 

A  comprehensive  theory  of  the  Cyclotron  Resonance  Maser  interaction  in  a 
waveguide  has  been  presented,  including  the  kinetic  theory  and  single  particle 
theory  approaches.  The  interaction  has  been  examined  for  both  TE  and  TM 
modes,  and  equations  have  been  expressed  in  a  simple  form  which  elucidates  the 
physics  and  facilitates  calculation.  The  kinetic  theory  was  used  to  derive 
dispersion  relationships  for  generalized  beam  parameters  and  to  calculate 
growth  rates  for  a  cold  beam.  The  single  particle  theory  was  used  to  obtain 
linear  and  nonlinear  results  for  a  cold  beam.  The  present  results  confirm 
that  the  form  of  the  equations  of  motion  and  the  dispersion  relationships  is 
the  same  fcr  TE  and  TM  modes  in  the  CARM  and  gyrotron  regimes.  The  theory 


provides  the  starting  point  for  treating  electron  beam  self-field  effects,  the 
use  of  plasma  filled  waveguides,  and  the  treatment  of  time  dependent 
phenomenea. 

The  calculations  carried  out  for  the  oscillator  configuration  as  a 

function  of  the  parameter  b  apply  to  essentially  all  possible  CRM 

configurations  including  gyrotrons  (b  ■  0),  CARM's  (b  •  0.5),  and  highly 

relativistic  doppler  shifted  configurations  with  high  pitch  angle  beams 

(b  ■  1).  Comparison  of  the  latter  two  regimes  shows  that  normalized 

efficiency  is  higher  at  b  -  0.5  than  at  b  ■  1 ,  but  the  single  particle 

efficiency  i3  higher  at  S  •  1  (n  „  £  100%)  than  at  b  -  0.5  (n  S  50%)  .  Thus 

sp  sp 

the  actual  efficiency  is  similar  in  the  two  cases  (20-25%). 

These  calculations  show  that  the  doppler-shifted  cyclotron  maser 
interaction  is  an  attractive  candidate  for  use  in  high  power,  high  frequency 
sources.  Nonlinear  efficiencies  of  order  25%  are  feasible  and  this  can  be 
increased  by  enhancement  techniques.  The  CARM  regime  appears  to  offer  an 
attractive  combination  high  interaction  efficiency,  large  doppler  upshift,  and 
low  beam  orbit  pitch  angle.  Examples  of  amplifier  and  oscillator 
configurations  are  given.  These  examples  show  the  potential  of  the  CARM 
configuration  as  a  high  power  millimeter-wave  and  sub-millimeter-wave  source 
operating  at  the  50-200  MW  level.  The  present  calculations  do  not  account  for 
beam  temperature  or  velocity  spread.  CARM  amplifier  simulation  studies  by  Lin 
and  Lin  (1985)  indicate  that  less  than  1%  axial  velocity  spread  is  required  to 
achieve  high  efficiency  in  a  configuration  similar  to  the  present  94  GHz 
example.  The  shorter  interaction  length  of  the  CARM  oscillator  should  render 
it  less  sensitive  to  velocity  spread:  the  calculations  of  Xanavets  and  Klimov 
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( 1 975 )  suggest  the  transverse  velocity  spread  (standard  deviation)  should  be 


less  than  5S. 
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Table  I:  100  GHz  CARM  Oscillator  Parameters 


* 

* 


Electron  Beam  Energy 
Optimum  Beam  Current 
Threshold  Current 
Initial  Transverse  Velocity 


600  kV 
250  Amp 
60  Amp 
0.46c 


Initial  Axial  Velocity 


0.76c 


Operating  Mode 


Cavity  Length 


4.6  cm 


Cavity  Wall  Radius 
Annular  Beam  Radius 


1 .2  cm 
1.1  cm 


Wave  Phase  Velocity 


1.13c 


Cavity  Output  Reflectivity 
Cavity  Q 


a 


Magnetic  Field 
Interaction  Efficiency 
Output  Power 


23  kG 


37  MW 


W! 


Table  II:  91*  GHz  CARM  Amplifier  Parameters 


Electron  Beam  Energy 
Optimum  Beam  Current 
Initial  Transverse  Velocity 
Initial  Axial  Velocity 
Operating  Mode 
Interaction  Length 
Waveguide  Wall  Radius 
Solid  Beam  Radius 
Wave  Phase  Velocity 
Magnetic  Field 
Saturated  Efficiency 
Output  Power 

Small  Signal  Bandwidth  (FWHM) 
Large  Signal  Bandwidth  (FWHM) 


1  MeV 
500  Amp 
0.33c 
0.76c 
TE11 
54.5  cm 
Q.5H  cm 
0.3  cm 
1 .015c 
13.2  kG 
22  % 
110  MW 
10  i 
3  * 


Table  III:  2  THz  Second  Harmonic  CARM  Amplifier  Parameters 


Electron  Beam  Energy 
Optimum  3eam  Current 
Initial  Transverse  Velocity 
Initial  Axial  Velocity 
Operating  Mode 
Interaction  Length 
Waveguide  Wall  Radius 
Annular  Beam  Radius 
Wave  Phase  Velocity 
Magnetic  Field 
Saturated  Efficiency 
Output  Power 

Large  Signal  Bandwidth  (FWHM) 
Average  gain  rata 


1.2  MeV 
1  kA 
0.3c 
0.9c 

TE10,1 
131*  cm 

0.5  cm 
0.4  CD 
1.0016c 
115  kG 
19  i 
223  MW 
0.04  % 

0.4  dB/cm 


• 
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Figure  5:  a)  Normalized  efficiency  for  first  harmonic  as  a  function  of 

normalized  current  for  b  -  0.454  and  FQ  -  0.0018.  The  results  have  been 
optimized  with  respect  to  &  and  interaction  length  p  .  b)  Peak 
normalized  wave  amplitude  corresponding  to  efficiency  shown  in  a),  c) 
Optimum  interaction  length  corresponding  to  efficiency  shown  in  a) . 


Figure  6:  a)  Normalized  efficiency  dependence  on  the  input  wave  amplitude 

A 

for  the  first  harmonic  (s«1).  The  normalized  current  is  I  -  0.60  and 

o 

D  •  0.47.  The  results  have  been  optimized  with  respect  to  A  and  u 


as  in  Figure  5.  The  optimum  values  of  saturation  wave  amplitude  and 
interaction  length  are  shown  in  b )  and  c). 


velocity  (and  frequency)  in  the  design  waveguide  (ru  =*  0.5*1  cm). 
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A 

Figure  8:  The  small  signal  (dashed  curve)  and  nonlinear  (solid  curve) 
first  harmonic  efficiencies  plotted  as  a  function  of  the  detuning 
parameter  A.  The  interaction  parameters  are  given  in  Table  II.  The 
small  signal  curve  has  been  multiplied  by  a  factor  of  10  for  clarity 


